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A  LINEAR  PROGRAMING  FORMULATION  OF  A  SPECIAL  QUADRATIC 
ASSIGNMENT  PROBLEM 

Abstract 

A  special  quadratic  assignment  problem  is  shown  to  be  equivalent 

3  2 

to  a  linear  programing  problem  with  n  constraints  and  n  variables  where 
n  is  the  number  of  elements  Co  be  assigned.  A  labeling  algorithm  similar 
to  that  for  the  linear  transportation  problem  ia  presented  for  solving 
the  problem.  An  example  is  presented  that  deals  with  "tr iangularizing" 
input-output  matrices. 
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Introduction 


This  paper  shows  that  a  special  quadratic  assignment  problem  can 
be  reformulated  as  a  linear  programming  problem. 

The  general  quadratic  assignment  problem  can  be  written  as 


max  T*  y'AY 

St  =»1J  = 

1 

^yU 

a 

L 

(1) 

0,1 

where  the  matrix  dot  product  is  defined  as  T»C  *  £  t  c  ,  i.e., 


<i,j> 


ij  ij 


is  sum  or  products  of  all  paired  indices,  see  Lawler  [3  ].  In  Section  1, 
we  show  that  for  a  special  form  of  T,  (1)  can  be  reformulated  as  a 
linear  programming  problem,  i.e.,  all  extreme  points  of  a  linear  set 
of  constraints  are  solutions  to  (1)  and  vice-versa.  In  Section  2,  « 

develop  a  labeling  algorithm  for  solving  the  linear  programming  formulation. 
In  Section  3  we  discuss  an  Application  of  this  formulation  to  the  triangu¬ 
lar  ization  of  input-output  matrices. 


Section  1:  Linear  Programming  Equivalent 

Consider  the  quadratic  assignment  problem  (1)  where  T  is  an  upper 
triangular  matrix  of  ones,  i.e., 


....  0  1 


T 


I  0 
\  o 


0  0 


This  quadratic  assignment  problem  can  be  viewed  as  finding  a  permuta¬ 
tion  P  -  (r,,...r  )  of  the  rows  and  columns  of  A  such  that  the  sum  of 
the  above  diagonal  elements  is  maximized.  That  is,  if  f^(P)  is  the 
value  of  (1)  for  permutation  P  Chen 


fA(P> 


n-1  n 

E  E  a  . 

i«l  j“i+l  rl,  j 


The  relationship  of  the  permutation  P  to  the  variables  Y  in  (1)  is 
that  Y  is  the  permutation  matrix  of  P,  that  is 


i  4  r 
i  *  r 


j 

j 


The  structure  of  the  problem  yields  the  following  easily  verified 
properties: 


Property  1:  For  any  permutation  either  or  a.^,  i  i  j 

must  appear  above  the  diagonal  so  that  either  or  a 

is  in  Che  function  f.  (P). 

A 


Property  2:  If  A  is  a  symmetric  matrix,  then  fA(P) “constant 


for  all  P. 


We  can  assume  by  Properties  2  and  3  that  a^  >  0  for  all  i  and  j. 
Let  us  consider  the  graph  G  ”  (X,U)  where  X  is  a  set  of  n  vertices  and 
U  is  the  set  of  directed  arcs  U  *  {(i»J)  i  4  jl.  We  will  associate 
with  each  arc  (i,J)  of  G  a  cost  a£j*  l*t  r  be  the  set  of  partial 
graphs  formed  with  arcs  of  G  such  that  (l.jJsG'cT  »  (j»i)  p  G* . 

Let  h(G’)  »  V  a  for  G'«r,  G*  »  (X',0).  If  we  consider  the 
(i.J)«C*  Li 

following  problem 

min  h(G') 

st  G'er  (2) 

and  G1  has  no  circuits, 

then  (1)  and  (2)  are  related  by  the  following  theorem. 

Theorem  1.  If  i  is  the  optimal  solution  te  (1)  and  w  the  optimal 
solution  to  (2),  then  z  ■  w 

Proof :  In  order  to  prove  this  theorem  we  will  show  that  every  feasible 

graph,  G',  in  (2)  corresponds  to  a  permutation  P'  in  (1)  and  that 
there  exists  no  permutation  in  (1)  for  which  this  correspondence  does 
not  hold. 

Lemma  1:  If  G'*r  then  (i,J)  f  G'  «*  (J,i)  e  G'  for  i  4  J. 

Proof:  Consider  the  sets  of  arcs  ■  {(i»j)>(J»i)  J  i  >  jl;  There 
are  n(n-l)/2  such  sets  and  by  hypothesis  only  onu  «rc  be  u.jvr. 
from  each  set  in  constructing  G'.  If  there  exists  a  set  such  chat 
S^nu'  »  0  then  |U'|<  n(n-l)/2  contradicting  the  face  that  G'cT. 
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Lemma  2:  If  there  exists  a  circuit  in  G ' cT  of  length  g  >  3  then  there 
exists  a  circuit  of  length  3. 

Proof:  Let  (i^.i^, . . . ,  i^^^.i^)  denote  the  circuit  G'eiT.  Now  suppose 
(ij,iq_2)f  G’  since  otherwise  we  are  done.  By  Lemma  1,  ( i^_2> ij) eG' • 
thus  there  exists  a  circuit  of  length  g*l.  This  lndirective  step  can 
be  repeated  to  prove  the  lemma. 

let  us  define  the  out  degree  N(i)  of  a  vertex  i  in  G'  as  N(i)=  £  1 

(i.j)eG' 

Lemma  3:  If  N(i)  -  N(j)  i  4  j  then  G'  has  a  circuit. 

Proof:  Let  N(i)  *  NCj^  »  K  and  let  (i, jj) ,(i, j2> . (i.j^eG'  . 

Now  suppose  (JjjtJeG'  and  t  4  t  «  2,...,K,  then  (t,i)eG'  by  Lemma  1. 
Therefore  t  «  i  *  2,...,K,  but  there  are  only  K-l  such  arcs  and 
since  N(jj)=K  there  Is  one  arc  ( J ^ , t)  for  which  t  4  j^,  i  *  2,...,K 
aud  there  is  a  circuit. 

Proof  of  Theorem  1. 

Let  G  be  feasible  in  (2).  Then  the  nodes  of  G1  can  be  ordered  such 
that  i  <  j  if  N(i)  >N(J).  Let  this  ordering  be  P'  «  (r,  ,r,, ,. . .  ,r  )  i.e., 
rj  “  k  if  N<*>  »  "  -  j.  Since  0  <  N(i)  <  n-l.N^)  »  N(i  )  +  1,  and  thus, 
P  is  a  peroutation  of  the  n  integers  (0,1, .. . ,n«l) ,  which  is  isomorphic  to 
the  peroutation  of  the  integers  (l,2,...,n).  Consequently  P'  is  a  permuta- 
tion  for  (1).  Notice  that  any  permutation  in  (1),  introduces  the  ordering 

on  G‘  by  the  out  degree  values.  To  show  that  z  -  w,  we  note  that  for  the 

r.  .  n 

permutation  ?*»*■£  Z  *  _  while  for  the  associated  graph  G'*.‘, 

i-1  j-i+1  rirj 

n-1  n 

w-h(G ')  -  T  a  -  v  r  a 

(i,J)eG’  i-1  J-l  rirj 

Q  •  £  •  D  • 
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Coroliarv:  The  optimal  solutions  to  (1)  and  (2)  are  related  as  follows 


ij 


1  if  N(i)  -  n-j 
0  if  otherwise 


Proof:  Obvious  from  the  proof  of  Theorem  1. 

Prom  Lemma  2  we  can  reformulate  (2)  as  followa: 

min  h(G’) 
st  e'er 

0‘  has  no  circuits  of  length  3. 

If  we  let  x^  =  1  if  <i,j)  e  G* 

=■  0  otherwise, 

then  the  following  aquations  do  not  allow  circuits  of  length  three: 

*o *  xJk  +  *ki  - 2  1  *  -* k 

V +  Ktk 

in  addition.  Property  1  stipulates  +  x^  *  1  i  4  j.  Thus  the 
quadratic  assignment  problem  has  been  reduced  to  the  following  integer 
linear  programming  problem 


n 

n 

£ 

i-l 

Tax 
j-1  « 

x .  . 

+  X  ..  +  X  .  < 

ij 

jk  ki  ~ 

X,. 

+  x.  ,  +  x  ,  < 

ik 

kj  ji  “ 

x  . 
1J 

4*  x  s  1 

*jk  * 

Xii 

*  0  for  all  i 

w 


i  *  J 


6. 


We  now  use  the  following  theorem  of  Bowman  1  1  ]  to  reduce  (4) 
and  correspondingly  the  quadratic  assignment  problem  (1)  to  a  linear 
program. 


Theorem  2:  The  extreme  points  of  the  polyhedron  formed  by  me  constraints 


x.  .  +  x  ..  +  x  <2 
13  jk  ki  - 


i  ¥  k  4  j  4  l 


Xik  +  \j  +  xji  -  2 


x .  +  x *  1 

i.l  jk 


1  ¥  J 


x .  .  >0 
J-J  " 


x .  .  =  0  for  all  i 
ii 


are  all  integer. 


The  result  of  this  theorem  is  that  wa  can  drop  the  integer  con¬ 
straints  to  (4)  which  reduces  the  problem  to  a  linear  programming 
problem.  However,  there  are  n(n-l)  variables  and  n(n-l)/2  +  n(n- 1) (n-2) /3 
constraints  which  causes  problems  with  solving  by  conventional  L.P. 
methods  as  n  grows  large.  In  the  next  section,  we  develop  a  labeling 
technique  to  solve  (4) . 


Section  2:  Algorithm 

In  this  section  we  describe  a  Labelling  type  algorithm  for  solving 

the  linear  programming  problem  (4).  This  algorithm  corresponds  to  a 

dual  simple::  routine,  but  the  amount  of  computer  storage  required  is 

smaller  than  the  standard  linear  programming  requirements. 

We  first  note  that  the  constraint  x  +  x  =  1  can  be  replaced 

ij  ,]i 

by  upper  bound  constraints  on  one  of  the  variables.  There.fore,  let 
us  reformulate  (4)  as 


u-1  n 

max  Z  Z  ( a  .  .  -  a  . . ) x  . 

i-1  >1+1  iJ  U 


x.  .  +  X..  -  X.,  •'  i 

ij  ik  ik  — 

i 

<  j  <  k 

(5) 

-x  .  .  -  x  +  x <  C 

ij  jk  ik  — 

i 

<  j  ■ :  k 

0  <  x .  .  i 

-  ij  - 

i 

'  j  • 

We  have  thus  reduced  the  problem  to 

one 

invo ivrng 

— variables, 

but  there  are  sLil.l  n(n-i)  (n-d)/3  constraints  excluding  the  upper 
bounds.  This,  of  course,  implies  solving  the  dual  of  (5).  Let  (5) 
have  the  matrix  representation 

max  cx 
s. t .  Bx  <  e 
and  -Bx  <  0 

°  <  x  <  1 

then  the  dual  to  (b)  is 

miu  ew  t  ev 
s.t.  B'w-B'u+  Iv> 
w,u,v,  >  0 


c 


■8- 


or 

min  ew  +  ev 

a.t.  -B'w  +  B'u  -  Iv  +  Is  ■  -c  (6) 

s,w,u,v,  >  0 

Now  since  the  columns  o£  B*  correspond  to  the  constraints  o£  (.5)  every 
column  can  be  denoted  by  a  triple  (i,j,k)  where  i  <  j  <  k.  Of  course, 
identifying  a  column  of  B'  is  equivalent  to  identifying  variables 
in  (6)  and  we  therefore,  use  the  following  notation: 


denotes  variable  u  .  i  <  j  <k 

lJK 


and 


(7a) 


denotes  variable  w^^  i  <  j  < 


In  a  similar  manner  we  let 


(0 , i , J)  denote  variable  s^  i  <  j 

and  (7b) 

(0,i, j)  denote  variable  v  i  <  j. 

Observe  that  the  variable  designation  also  allows  one  to  generate  the 
appropriate  column  representation  in  (6).  Because  of  this  fact  we 
will  develop  an  algorithm  that  keeps  track  of  the  basis  elements  of 
(6),  but  not  their  matrix  inverse.  Then  when  we  decide  to  enter  a 
new  variable,  we  will  find  its  representation  by  a  labelling  algorithm. 
This  is  quite  similar  to  the  approach  used  in  transportation  algorithms. 

In  order  to  determine  Incoming  variables  (or  optimality)  we  will 
use  the  fact  that  if  we  know  the  dual  solution  associated  with  a 
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basis  of  (6),  i.e.,  values  for  x.^,  Chen  the  reduced  costs  for  w  and 

u  are  the  slacks  of  the  respective  constraints  for  x,  .  +  x.,  -  x.t  <  1 

ij  jk  ik  - 

and  -x.  .  •  x,.  +  s„  <0  and  the  reduced  costs  for  v  and  s  are  re  spec - 
ij  jk  ik  - 

tively  1-x^  an^  Xij*  ^en  these  reduced  costs  are  all  positive  we 
are  optimal,  otherwise  we  introduce  the  appropriate  variable  into  the 
basis.  In  order  to  understand  the  following  algorithm  we  again  point 
out  the  variable  designations  (7a)  and  (7b)  allow  for  the  generation 
of  the  appropriate  column  of  (6)  and  thus,  there  is  no  need  to  ever 
have  the  complete  representation  of  (6). 

We  denote  as  follows  various  storage  units  used  in  the  algorithm. 
X:  Upper  triangular  matrix  whose  values  x_  correspond  to  the 
solution  of  (5)  associated  with  a  basis  of  (6)  i.e.,  they 
are  the  reduced  costs  associated  with  variable  s. 

C:  Right-hand  side  values  associated  with  current  basis  of  (6). 
B:  Upper  triangular  matrix  b  is  the  identification  of  the 
basic  variable  associated  with  row  (i, j)  of  (6) . 

S:  n(n-l)/2  lists  denoted  S^.  List  contains  the  basic 

variables  that  have  non-zero  entries  in  row  (i,J)  of  (6). 

The  other  elements  used  are  appropriate  temporary  storage  units. 
Step  :  Initialization 

C1J  ■  Kj  •  a)ib 

li  (aLj  -  ajl)  <  0,  x.j  -  0,  btj  =  (0,i,j),  S.j  =  l(0,i,J)} 
Otherwise,  *  1,  b  -  (0,i,J),  S_  =  {(0,t,J)} 

Remark:  It  is  easy  to  confirm  that  this  produces  a  feasible  basic 

solution  to  (6)  made  up  of  the  variables  v  and  s. 
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Optimality  Check 


If  0  <  x.  .  <  1  for  all  i  <  j  and 

-Lj“ 

0  5  ^  +  xjk  •  x,k  <  1  Lor  all.  1  <■  j  <  k, 
stop,  this  solution  is  optimal. 


then 


Remark:  Let  v  =  {(r,s,t)  or  (r,s,t)j  he  the  entering  variable  corresponding 
to  violation  of  one  of  the  above  restrictions.  Note  r  1  0  is 
permissable.  Go  to  Step  2. 

Remark:  When  r  =  0,  addition  on  index  0  implies  no  operation  and  is 
used  only  for  notational  convenience. 


Step  2: 


2A:  Let  z..  =  |S.  .|  (the  number  of  elements  in  list  S.  .)  . 

tj  ij  i-J 

Let  z  -  z  'tl,  z  —  z  d*l,  z  ~z  d*l. 
rs  rs  rt  rt  st  st 


2b:  Reduction  Step 

2B1:  D  =  0  (D  is  a  list) 

2B2 :  If  ?. ,  .  =  1  tor  any  z .  . ,  let  P  -  SJ  .  -  S,.  0  D  » 

ij  7  J-J  i3  ij 

-  [(a,b,c)  or  (a,b,c)}  and  go  to  Step  2B3 . 

Otherwise  go  to  Step  3. 

2B3:  D  -  D  U  P  and  z  ,  =  z  ,  -  1,  z  *  z  -1  and  z,  =  z,  -  1. 

ab  ab  ao  an  be  be 

Go  to  2P2. 

Remark:  Step  2B  eliminates  basic  variables  that  cannot  be  used  to 

represent  (r.s,t)  or  (r,s,t)  because  of  having  only  one  non¬ 
zero  clement  in  a  row  of  the  basis  matrix  that  has  a  zero  in  the 
entering  column. 


1. 


T  I  Jj'  ~  I 

3N:  II  v  -  (r,s,t),  '-rs  ’  Jrt 

Step  3C . 

Otherwise .  if  Krs  " 

Co  to  3C . 

3C :  Let  Q  -  Un  ordered  list)  . 

3D: 

T  -  !  1-id  E  4  o  tor  SO! 

3D1:  If  ,  =  '  1  iUld  Lii  T 


-1.  Go  to 


■U  E  -  F 


..  -  i  and  E.  .  4  Q  for  some  i  -  j»  Let  ?  '  Sij  Sij 

*  i  j  1 1 


Let  b  =  r  and  L. (b  )  -  +1* 
pq  -  1  H 


i.  to  3E. 

3D2: 

It  E.. 

=  0 

tor 

alL 

i.r> 

feo 

to  Step  4 

• 

3D  3: 

If  *ii 

=  G 

tor- 

some 

Eij 

4 

0.  Go  to 

Step 

3 . 

it  zL. 

L  ?. 

tor 

all 

U.  .  1 

11 

l  0 

II 

fc- 

some 

element  of 

3 .  .  " 
iJ 

3 .  . 

i.1 

n  d 

tor  E 

ij  4 

0. 

Let.  b 

pq 

=••  V. 

L.(b  )  - 

i  pq 

Go  to 

3E . 

It  b  ~  (a,b,c)  then  either  (a,b)  =  (i»j)  and  s X 

L2(bpq)  =  sign  (-E.  .)  or  (a,c)  •  d»i>  and  set 

L  (b  )  --  sign  (E..)  or  (b,c)  ^  (i,j)  and  set 
2  pq 

L0fb  )  -  sign  (-E .  .)  ? 

2  pq 


E  ,  -  E  .  +  E  (b  )  •  5 
ab  ab  2  pq 

3  -  V  -  E„(b  )  '  1 

ac  ac  2  pq 


p  -  E  +  L  <h  )  '  L  r'°  ,0  3F‘ 
l'hc  be  2  pq 

.  •  -  c  k  -  (a  b  e) then  either  (a>b)  -  (i»)) 

-rri*  otherwise  if  b  -  la.o.wwcu 

—  PS 

and  set  L,^)  *  sign  <K5j>  0r  'a,c)  =  (“’J>  ^  *** 

L  (b  )  =  sign  (-E  .)  or  (b,c)  =  d.j)  and  set 

?  pq  JJ 

4(bpq)  '  9lRn  (Eij>  '• 


ab 

Kab  " 

S(\q> 

dC 

E  ^ 
nr 

W 

Ev  •' 
be 

F'bC  " 

W 

Go  to  3E 
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3F :  z  =2  -  1,  z  =  z  -  1,  zt  .  =  zt  -  1,  D  =  D  0  P 

ab  ab  ac  ac  be  be 

Q  -•  Q  U  (b  ,  L.(b  ),  L  (b  ))  Go  to  3D. 

pq  i  pq  2  pq 

Remark:  Step  3  finds  the  representation  of  the  incoming  column  v 

in  terms  of  the  present  basis.  is  a  label  that  is  one 

if  there  is  only  one  non-zero  entry  in  a  row  which  has  a 
non-zero  value  in  v  or  its  subsequent  representation. 

Step  3t)l  always  tries  to  force  a  permanent:  labelling  since 
this  is  a  forced  relationship, 

I,  is  a  label  that  gives  the  sign  of  the  basis  element 
b  in  its  representation  of  v.  Since  the  set  is  unimodular, 

pq 

we  know  the  coefficient  of  a  basis  vector  is  0  or  +  1, 

Finally,  we  note  that  3P2  stops  our  labelling  procedure 
when  wt-  have  found  a  linear  combination  of  the  basis  vectors 
that  yields  v. 

Step  4:  Change  of  Basis. 

4A  :  Let  Q+  he  the  set  of  Q  such  that  L„(b. .)  =  +  for  b. ,6  Q  and  0 

2  ij  13 

be  the  set  of  Q  such  that  L  (b..)  =  -  for  b..£  Q .  Let  min  c..  - 


2  ij 


c 

and 

w  = 

b  . 

pq 

pq 

c .  . 

-  c 

.  + 

c 

for 

all 

b.  .  ' 

i.l 

IJ 

pq 

rj 

c  .  . 

-  c 

c 

for 

all 

b..  1 

ij 

ij 

pq 

ij 

b 

=  V 

and  S 

=  S 

Uv, 

S  _ 

pq 

r 

a  r 

s 

rt 

and 

if 

W  = 

(a,b 

,c)  or 

■  (a, 

b  ,c) 

S 

_ 

S 

-  W  and  S. 

=  S. 

ac 

ac 

be 

be 

For 

(0, 

i.j) 

6  B, 

set  x 

•  e 

ij 

o. 

b.  ar* 
ij 


ab  ab 


Co  to  4C . 
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4C:  For  (i,j,k)  €  B,  i  4  0  set,  =  0  and  for 

(i,j,k)  (r  B  set  +  x  -  x^  ~  1»  Solve  these  constraints 
with  the  conditions  of  4B  to  obtain  new  X  matrix.  Go  to  Step  1. 
Step  5:  Back-tracking. 

5A:  Back-track  in  list  Q  to  the  last  element  labelled  0  (temporary) 

say  b  .  All  elements  encountered  in  the  back-tracking  prior 
mn 

to  b  say  b  must  have  L, (b  )  =  +1  and  are  treed  i.e.,  for 
mn  pq  1  pq 

all  p,q.  Perform  5AI,  5A2,  and  5A3. 


541:  If  b  =  (a ,b ,c) ,  then 

pq 


E  =  E 
ab  ab 


L,,(b  )  *  1 
/  pq 


E.,c*  Eac+  L2(V 


Ebc  *  Ebc 


L  (b  )  '  1.  Go  to  Step  5A3 . 
2  pq 


5A2:  Otherwise,  if  «  (a,b,cj,  then 

E„b  -  Eab  +  h(W  •  1 
Eac  ■  Eac  -  VW  ’  1 


‘be 


=  Ebc+  Go  t0  Step  5A3. 


2'  pq 


5A3 : 


D  =  D-b,Q=Q  -  (b  ,L  (b  ),  L  (b  )) 
pq  pq  1  pq  2  pq 

2,  =  Z  ,  +  1  .  2  =2  +  1,  Z.  =  Z,  +  1. 

ab  ao  ac  ac  be  be 


5B :  L. (b  )  =  +1.  Go  to  3D. 

1  mn 

Remark:  Note  that  the  only  change  for  b  is  that  it  is  permanently 

labelled  rather  then  temporary.  This  makes  sure  that  it  will  be 
freed  when  we  back-track  to  some  point  prior  to  Its  assignment. 


While  this  algorithm  may  seam  quite  involved  compared  to  simple 
simplex  operations,  one  should  note  that  its  advantage  comes  in  its 
small  core  requirement  compared  to  a  linear  programming  approach. 

If  we  assume  that  the  original  matrix  is  generated  rather  than  stored, 


our  core  requirement  for  simplex  solution  of  (b)  is  essentially  that  of 

2 

the  basis  inverse.  This  is  or  on  t^e  or^er  Df  n^.  For  the 

above  algorithm,  one  needs  to  store  elements  for  X,C,b,Z,E,  and 

5  /  «  K 

2  labels  and  a  maximum  of  for  S,  that  is  approximately  5n(n-l). 

For  large  n  this  savings  is  considerable.  For  n  *  20  the  L.P.  may  re¬ 
quire  over  36,000  storage  spaces  while  the  proposed  algorithm  requires 
roughly  1900. 

In  the  next  section,  we  discuss  the  application  of  this  algorithm 
to  a  quadratic  assignment  problem  involving  the  triangular izac ion  of 
Input-Output  matrices. 
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Section  3:  Example  -  Triangular irat Ion  Input-Output  Matrices 

The  permutation  of  rows  and  columns  of  an  input-output  matrix, 
denoted  A,  that  maximize  the  sum  of  the  above  diagonal  elements  provide 
a  hierarchical  listing  of  the  sectors  of  the  economy  as  "users"  and 
"suppliers". 

For  the  economic  interpretations  and  a  heuristic  method  for  triang 
ularization,  see  Simpson  and  Tsukui  {6J.  For  an  implicit  enumeration 
algorithm,  see  15). 

It  is  easily  verified  that  the  problem  of  permuting  the  rows  and 
columns  of  A  makes  it  a  quadratic  programming  problem  and  since  we  are 
concerned  with  the  sum  of  the  above  diagonal  elements  we  have  the  T 
matrix  as  in  Section  1. 

The  following  example  is  a  5x5  subset  of  an  original  37x37  input- 
output  matrix  presented  by  Leontief  [4  }. 


Agrlc.  and  Fish 

2453 

3896 

2195 

15 

317 

Food  and  Kindred  Prod. 

538 

1427 

61 

8 

0 

Textile  Mill 

14 

0 

1321 

2913 

0 

Appatel 

9 

50 

0 

1471 

0 

Lumber  and  Food 

34 

25 

20 

0 

1817 

16 


1 


f  (0,1,2)  (0,1,3)  (0,1,4)  (0,1*5)  y 

(0,2,3)  (0,2,4)  (0,2,5) 

‘  (0.3,4)  (0,3,5)  j 

V  (0,4,5)  / 


S  =  ». 

Step  1;  x23  +  x ^  -  x24  »  2  and  v  »  (2,3,4) 


Step  2: 


2A : 


X  » 


2Bl: 

2B2: 

2»3: 


2B2: 

2B3: 


D  *  P 

z1>2  »  1,  P  -  (0^2) 
o  *  1(oTU2)} 

/0  1  1  1\ 


*1,3  *  l»  P  *  <°*l*3> 

D  -  ((W).  (0X3)} 


2B2: 

2B3: 


(0,1,4) 


0  -  l (0,1, 2),  (0,1,3), 


(0,1,4) 


f  0  0 

2 

\ 


0 

2 

2 


1 

1 

1 

i 


\ 


I 

I 


/ 
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2B2 

and  2B3  repeat  for  5,z->  5**3  5*  and 

,  ending  with 

2B3: 

:  0  =  l (0,1 ,2) ,(0,1 ,3) , (0,1,4) , (0,1,5) , (0 

,2, 5),  (0,3, 5),  (0,4, »} 

000 

0 

2  2 

2  *  1 

0 

2 

0 

0 

0  0  0 

0  0  0  0 

0 

3A: 

E  =  00 

0  11 
z  ■ 

0 

0 

0  1 

0 

V 

0 

0 

0  0  0 

0 

3B: 

E  *  1  -1 

0 

i 

0 

0 

3C: 

Q  =  * 

3DI 

'*  *2.3“  1  and  E2 

y  0  p  -  (57575). 

Ll(b2,3)  *  +1 

3E: 

L2(b2>3)  =  + 

Q  *  U*>2,3’  1>  +)5 

000 

0  0  0  0 

0 

E  -  0  -l 

0  0  1 
z  ■■ 

0 

1 

0  1 

0 

0 

0  ' 

3D1 

:  Zj  ^  »  1  and  E 

2>4  *  0  p  -  (0,2,4) 

3E: 

h(h2.i  ‘  +•  Q 

*  l(b2  3, !>.•+)»  (b2 

0  0  0 

0  0  0  0 

0 

E-  00 

0  0  0 

z  ■ 

0 

1 

0  l 

0 

0 

0 

18 


19 


The  optimal  solution  has 

/Kill 

j  o  X  1  1 

**  o  0  x  1 

|  0  0  0  X 

\0  1  1  1 


i  1  *  2 


I 

i  !-  l 

i 

1  !  o 


sc  that  the  optimal  arrangement  is  (1,5, 2, 3, 4)  i.e., 


Agric.  &  Fish 

/2453 

i 

317 

3896 

2195 

15 

Lumber  &  Food 

34 

1817 

25 

20 

0 

Food  &  Kindred  Prod. 

538 

0 

1427 

61 

8 

Textile  Mills 

14 

0 

0 

1321 

2913 

Apparel 

\  9 

0 

50 

0 

1471 

or  in  reduced  form 


X 

283 

3358 

2181 

6 

0 

X 

25 

20 

0 

0 

0 

X 

61 

0 

0 

0 

0 

X 

2913 

0 

0 

42 

0 

X 

While  the  algorithm  may  seem  quite  cumbersome  for  such  a  small 


problem,  efficient  use  of  hit  processing  routines  and  proper  programming 
techniques  greatly  simplify  the  updating  procedures. 
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Couc lusion 

We  have  shown  that  a  special  form  of  the  quadratic  assignment 
problem  can  be  reformulated  as  a  linear  programming  problem,  and  that 
because  of  the  structure  of  the  linear  program  an  algorithm  can  be 
developed  that  parallels  the  simplex  method  but  at  a  greatly  reduced 
storage  requirement.  This  might  indicate  that  other  quadratic  assign¬ 
ment  problems  may  have  equivalent  linear  programming  formulations  lor 
specific  T  matrices  in  the  formulation  (1) .  These  of  course,  will  noc 
necessarily  have  to  correspond  to  integer  solutions  as  in  the  case  or 
the  in  the  above  problem,  but  rather,  may  Call  in  the  category  of 
searching  for  appropriate  permutation  polyhedra  which  are  discussed  by 
Bowman  [1],  One  should  further  note  that  the  algorithm  of  Section  2 
is  really  dependent  on  the  constraint  set  of  (5)  and  not  its  objective 
function  and  th\is,  is  applicable  to  any  problems  where  this  constraint 
set  arises  whether  or  not  it  comes  from  the  Quadratic  Assignment  problem 
(1)  .  For  a  discussion  of  this  constraint  set  in  the  context  of  transi¬ 
tivity,  see  Bowman  and  Coiantoni  (2j. 
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